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Abstract. The overall solid-to-solid phase transformation kinetics under non- 
isothermal conditions has been modelled by means of a differential equation 
method. The method requires provisions for expressions of the fraction of the 
transformed phase in equilibrium condition and the relaxation time for transition 
as functions of temperature. The thermal history is an input to the model. 
We have used the method to calculate the time/temperature variation of the 
volume fraction of the favoured phase in the a f) transition in a zirconium 
alloy under heating and cooling, in agreement with experimental results. We 
also present a formulation that accounts for both additive and non-additive phase 
transformation processes. Moreover, a method based on the concept of path 
integral, which considers all the possible paths in thermal histories to reach the 
final state, is suggested. 



1. Introduction 

The kinetics of phase transformations in solids often involves the effects of heating 
and/or cooling rates. This is because first order phase transformations generally occur 
by concurrent nucleation and growth of the new phase and that both these mechanisms 
are time and temperature dependent [lj [2] . The overall phase transformation kinetics 
(nucleation plus growth) is represented by the fraction of transformed material in 
the system as a function of time and temperature 4> = (f>(t,T). This attribute 
under isothermal conditions will results in a timc-tcmperature-transformation (TTT) 
diagram, which is a complement to the phase diagram in equilibrium thermodynamics. 
Besides the nucleation rate and the growth rate, other factors that determine 4> include 
the density and distribution of nucleation sites, the overlap of diffusion fields from 
neighbouring transformed volumes, and the collision by adjacent transformed domains 
®. 

A phenomenological stochastic kinetic model for the overall phase transformation 
under isothermal conditions was solved exactly by Kolmogorov [4] and later 
independently by Johnson and Mehl [5] and Avrami [6], the KJMA model. Early 
treatments of non-isothermal transformation kinetics are found in papers by Avrami 
[7] and Cahn [5]. Avrami [7] showed that for a special case where the nucleation 
rate is proportional to the the growth rate of the favoured phase, over a temperature 
range, the non-isothermal transformations can be considered as a series of isothermal 
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reactions at every time step that can be linearly superposed to give the non-isothermal 
condition (additivity rule). Later, Cahn [5] argued that phase transformations that 
involve heterogeneous nucleation quite often obey an additivity rule. He noted that for 
such situations, a non-isothermal transformation can be related to an isothermal one 
through simple rate rules. More specifically, Cahn showed that if the transformation 
rate (d(f>/dt) depends only on <f> and temperature T, i.e. only on the state variables 
and not on the temperature path by which it had arrived to that state, then can be 
considered as an additive quantity. This situation occurs if the nucleation sites were 
consumed early in the reaction, i.e. site saturation occurs, and if the growth rate is a 
function of instantaneous temperature only. 

Cahn's additivity principle has been used by many workers for the evaluation 
of non-isothermal transformations in materials; e.g., austenite to ferrite+pearlite 
transformations in steels [9] I10| 111] and a similar kind of transformation in titanium 
alloys under cooling [12] H3J- A number of investigators have evaluated the 
applicability of the additivity rule in detail [21 [TSl EH H3 [TH] . A general modular 
approach, which accounts for the three phase transformation mechanisms, nucleation, 
growth and impingement applicable to both isothermal and isochronous conditions, 
have been developed, discussed and applied to phase transitions in alloys [19] [20] [21] . 
A recent overview of numerical and analytical methods for the determination of the 
kinetic parameters of a phase transformation is provided in |22j . An analytical solution 
for non-isothermal KJMA rate equation comprising separate activation energies for 
nucleation and growth when the transformation occurs under continuous heating has 
been obtained (23] . In an integral concurrent approach, Elder et al. [24] have studied 
the non-isothermal glassy or amorphous metals (such as Fe-B, Cu-Zr and Mg-Zn alloys) 
by solving a set of coupled physically-based equations for time-dependent spatially 
uniform temperature field, for the nucleation and growth of single crystallite and the 
Kolmogorov [4] formula for (f>, simultaneously. 

In this paper, we employ a method for calculation of the volume fraction of the 
new phase as a function of time and temperature during phase transformation in non- 
isothermal conditions. The method satisfies Cahn's additivity rule and assumes that 
the system is not too far away from equilibrium. It requires the specification of the 
fraction of transformed phase in equilibrium as a function of temperature. The method 
has been used to compute the phase transformation behaviour of a zirconium alloy 
under both slow and rapid heating and cooling (up to ±100 Ks _1 ). It is equivalent 
to the model for grain boundary nucleation when the nucleation rate is high and site 
saturation occurs early during the reaction [25| . The applicability of the method could 
be computations of transformation behaviour during fabrications subject to a variety 
of heat treatments, for example [26] , or in-service material performance under extreme 
conditions [27]. We shall also outline a general method for treating cases for which 
the volume fraction of the transformed phase is thermal history dependent. 

The organization of this paper is as follows. In section [21 a kinetic model for 
non-isothermal transformation, in differential form and in integro-differential form, 
is presented. The application of the differential method to the solid state phase 
transformation in a zirconium alloy is presented in section [3] In section [4] we 
discuss the appropriateness of the model both from a theoretical stance and empirical 
applicability. In section [5] we summarize the main results and remark on possible 
future directions. 
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2. Kinetic model 

One common approach to model the kinetics of non-isothermal phase transformation 
is to utilize the so called additivity rule [HE]. The rule may be stated as follows: The 
temperature history (temperature vs. time) is subdivided into a number of small type 
steps. Then the time spent for a volume element of the material in time interval At; 
at a given temperature Tj divided by the incubation time t X i at which the reaction 
has attained a certain fraction x of completion, can represent the fraction of the 
total nucleation time required for formation of the new phase. When the sum of 
such fractions reaches unity, the transformation starts to occur. Symbolically, it is 
expressed as 

« f *U = 1. (1) 
t[ txi Jo t*[T(s)] 1 ' 

A suitable material parameter for tracing the progress of solid-state phase 
transformation is the transformed volume fraction y as a function of time t and 
temperature T with the property < y < 1. 

Suppose the transformation is additive and the rate of transformation is only a 
function of the amount of transformation and temperature [8], namely 

^ = f(y,T), (2) 

The time it takes for a certain fraction y x of the new phase to be completed is 

t x = I"* f(y,T)- l dy, (3) 



Equation @ is a sufficient condition for the additivity, since it satisfies the additivity 
rule (HJ) through equation (J3]). 

Following |10j , we consider that y is not too far from its steady-state or equilibrium 
value y s (T) at a given temperature T, with f(y s ,T) = 0. Thus, the first term in a 
series expansion of / about y s yields 

dy y s (T)-y 

dt Tc (T) ' { ' 

where r c is a characteristic time of phase transformation and formally corresponds 
to t c = (df /dy\y— y ) . Note that equation (H} is controlled by two external 
temperature-dependent functions, i.e., t c (T) and y s (T), with y s {T) being the fraction 
of a new phase reached at temperature T after infinitely long time and < y s (T) < 1. 
Both these functions are material specific temperature-dependent quantities and 
can be deduced from experimental data on properties of a particular material or 
derived from appropriate models verified with such data. By putting (j> = y/y s and 
k(T) = t~ (T), we rewrite equation ((4]) in a reduced form 

^=fc(T)(l-0). (5) 



Integrating this equation gives 



(j){t) = 1 - exp - / k[T(s)}ds). (6) 



Under non-isothermal conditions (heating/cooling), (f> not only depends on t and 
T directly, but also on the way it has reached that temperature at a given point in 
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Figure 1. Time (t) temperature (T) histories. A sample undergoing a phase 
transformation when subjected to a thermal history from state (to, To) to (t n ,T n ) 
via different paths would reach different stages of transformations at (t n ,T„). 



time, see figure [TJ This process may formally be described by means of a generalized 
Langcvin-type equation. On these conditions, the temperature follows a path which 
may be an arbitrary function of time, and hence equation (|4|) needs to be solved 
numerically (section [3. 2[) . 

We next present a more general theoretical treatment of the kinetics of phase 
transformation. As pointed out in [8] , if the transformation rate depends only on y 
and T, i.e. on the state of the system rather on the thermal path by which it has 
reached that state, then the transformation is additive in the sense of equation (JTJ) . 
By making the ansatz /(</>, T) = h{T)g{4>), we write equation p} in the form 

^ = KT)g{cp). (7) 

As an example, for the KJMA model d Appendix A[ ), we write 

h{T) = TOfc(T) 1 " 7 (8) 

g{4>) =[-ln(l- (9) 

with 7 = (m — l)/m and m is the overall growth exponent \ Appendix A[ ). In our 
special case, m = 1, thereby (JT]) yields equation ((4]). 

We generalize the differential equation ((TJ to overcome the restriction of additivity 
[28] . In the manner of [29] , we write 

^=f M t T (*) - nmmw + m 

where M(t) is a memory kernel, which accounts for the effect of non-additivity and 
r](t) is a random force, e.g. a Gaussian white noise, which stands for the thermal 
fluctuations of the system during phase transformation. An important property of 
rj{t) is that it vanishes on the average, i.e. {rj{t)) = 0, and it is un-correlated to 4>{t), 
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i.e. (4>(t)r](t)) = 0, see e.g. [3D]. Also, M(t) is proportional to the spectrum of the 
random force, viz. 

M(t-t')^-^-(n(t)n(t')). (11) 

If now M(t) decays in a time tm and the applicable time is much longer than tm, 
then 

M(T[t] - T[t)]) = h(T)S(t - i), (12) 

where S(t) is the Dirac delta function. Thus, with relation (fT2"|) . Eq. (JTDJ) reduces to 
(J7J, that is when only the additivity rule is applicable. 

Let us now consider the case of non-isothermal transformation with a path- 
dependent characteristic. Suppose the volume fraction transformed obeys a general 
relation in the form 

^ = F[P(T)], (13) 

where Dy/DT is the functional derivative of y with respect to T and F[cp(T)] is 
a functional of tp in the sense that not only it depends on particular values of T 
and t, but on the function tp and all of T and t which are covered by tp, moreover 
< F[ip(T)] < 1. The form of this functional defines the transformation rule. The 
function tp may be described as a time integral of a rate of transformation k(T) [8], 131] 

tp[T(t)\ = f k(T(s))ds. (14) 



In case of interest to evaluate all the possible paths in the phase transformation 
from point (To, to) to (T, t) in the temperature-time plane, we re-express equation (|13|) 
in terms of a path-integral in the form 

f T,t 

y(T,t;T ,t ) = F[T(s),T(s)]d[T(s)}, (15) 

JT ,t„ 

where T = dT/dt and d[T(i)] denotes the product of infinitesimal steps in the thermal 
history, i.e., d[T(f)] = ^Tj, with dTj being the temperature increment at time tj. 

Hence, in the path- integral formulation, the time interval (to,i) is divided into n + 1 
intervals of equal length separated by time points t\, t^, ■ ■ ■ t n , t at which a material 
volume element (particle) is at temperatures T\, T 2 , . . . T n , T, respectively. Assigning 
a functional form for F by using a suitable model for isothermal phase transformation, 
a path-dependent description for the evolution of y under non-isothermal conditions 
is obtained. 

We choose a functional form for F[T(t),T(t)] in equation (fT5]) as 

F[T(t),t(t)]=6[T(t)]-exp{ - [ f k[T(s),T(s)]ds] m \ (16) 

where S[T(t)] is a functional Dirac delta distribution. Substituting equation p6[) into 
(|15p . with normalization, gives 

4>(T,t) = l- exp{- / k[T(s),T(s)]ds \d[T(t)}. (17) 

JTo,to J to 

This equation represents a path-integral description (or propagator) of the KJMA 
relation for the new phase development in non-isothermal conditions; and the 
functional integrand (|16[) may be interpreted as the probability-density for the new 
phase to follow a specific trajectory T(t), see figure Q] 
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3. Application 

3.1. Experimental data on zirconium alloys 

In this subsection, we make a short survey of experiments reported in literature, which 
we have used to select the input model parameters for the application of our model 
to zirconium base alloys. These experiments also provide data for model rctrodictions 
and predictions. We consider the kinetics of phase transformation of Zircaloy-4 (Zr- 
1.5Sn-0.2Fe-0.1Cr-0.12O, by wt%). Zirconium in solid state undergoes an allotropic 
transformation from the low temperature hexagonal closed-packed (hep) a-phase to 
body-centered cubic (bec) /3-phase at 1138 K [32] ■ On cooling, the transformation is 
either bainitic or martensitic depending on the cooling rate, with a strong epitaxy of 
the a-platelets in the former (3 grains [32] , 

Solid state phase equilibria of Zircaloy-4 have been investigated experimentally 
[33] . who reported a prevalence of four phase domains: a + \ up to 1081 K, a + (3 + \ 
from 1081 to 1118 K, a + [3 between 1118 and 1281 K, and /3-phase above 1118 K. 
Here, \ refers to the intermetallic hexagonal Laves phase Zr(Fe,Cr)2, see e.g. |34j . 
Quenching Zircaloy from /3-phase in moderate cooling rates produces two variants of 
Widmanstatten structure, namely, the basketweave and the parallel-plate structure 
[35] [55]. However, at cooling rates greater than 1000 Ks _1 a martensite structure is 
observed, while for very slow cooling rates, < 0.5 Ks _1 , the needle-shaped structure 
is rarely seen |36j . The overall transition in Zircaloy-4 has been studied by a 

number of workers in the past [33 HO [39l [13 EH [42] and more recently [27l [43] , which 
include also experiments on Zr-Nb alloys. 

Forgeron et al. [27] studied a <-» (3 transition of the Zr alloys by determining 
both their equilibrium (steady-state) temperature-dependence and their transient 
behaviour, with respect to the fraction of volume transformed, for the heating/cooling 
rates from ±0.1 to ±100 Ks -1 . More specifically, Forgeron et al. [27] determined the 
equilibrium behaviour of a//3-phase fraction as a function of temperature by means of 
calorimetry measurements. They used slow heating/cooling rates from 0.1 K/minutc 
to 20 K/minute. Moreover, they carried out direct measurements of the a//3-phase 
fraction by employing image analysis techniques on samples annealed for a few hours at 
different temperatures then quenched to room temperature. For kinetic measurements, 
they used dilatometry, where thermal cycles were applied on tubular samples, 12 mm 
in length, in vacuum or helium gas. As for calorimetric measurements, the uncertainty 
of the relative phase fraction measurement was less than 5%. 

3. 2. Computations 

In the model for computation of the relative phase fraction as a function of time and 
temperature, two functions y s (T) and r c (T) appear in equation j4|, which need to be 
specified. Let us consider first the former. The experimental data for the temperature 
dependence of steady-state volume fraction under phase transition suggest that y s {T) 
has an S-shaped or sigmoid form. For this reason, we have selected for the a ^> [3 
transition in Zr alloys a relation of the form 



for the equilibrium /3-phasc volume fraction at temperature T. Here, T cent and 
T S pan are material specific parameters that are related to the center and the span 




(18) 
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Figure 2. Temperature dependence of the characteristic time for Zircaloy-4 in 
the range of interest (top left) and its time-dependence for a triangular-shaped 
temperature pulse (bottom panel). 



of the mixed-phase temperature region, respectively. They are determined from the 
measured phase boundary temperatures T a and T@ through 



T n + T, 



0. 



Tp — T r , 



1079 K and 



1159 and T spa „ = 44 



2 , J- span — 23 ' V 1J ) 

Here, T a and Tp are defined as the temperatures that correspond to 99% a- and (3- 
phase fractions, respectively. For Zircaloy-4, the data in [44] gives T a 
Tp = 1273 K. We have used the data in [22|43] to obtain T cent 
K. 

The relation for r c (T) or its inverse the rate parameter, in general, k = r" 1 
depends on both the nucleation rate and growth rate of the new phase, which are 
strongly temperature dependent. It is usual to adopt an Arrhcnius-type relation for 
the rate parameter [3TJ [18] 

E 



k(T) = fco exp 



(20) 



k B T{t). 

where ko is a kinetic prefactor, E the overall effective activation energy and ks the 
Boltzmann constant. For Zircaloy-4, using the data on volume fraction [27], we 
found best fit values: k a = 60457+ 18129|Q| (s" 1 ) and E/k B = 16650 (K), where 
Q = dT/dt is the heat rate (Ks^ 1 ) in the range 0.1 < \Q\ < 100 Ks" 1 . Typical 
variations of t c (T) = k(T)^ 1 for a triangular temperature pulse is shown in figure [21 
Figurc[3]illustratcs the temperature rate dependence of r c as a function of temperature 
according to equation ([2"0f . 

The experimental results on Zircaloy-4 indicate that the starting temperature 
for the onset of phase transformation is temperature rate dependent [27] • This 



Transformation kinetics of alloys 



8 




Q( > LJ i ' 1 ' 

1100 1150 1200 1250 1300 1350 1400 
Temperature (K) 



Figure 3. Characteristic time as a function of temperature and temperature rate 
for Zircaloy-4 according to equation (50) . 

phenomenon has also been observed in other systems, e.g., fcrritic-pearlitic 
transformation in steels [1Q[ 111] and a <-> a + (3 transition of the titanium alloys 
[12]. In our modeling of this effect, for Zircaloy-4, we have related the onset of 
a — > a + (3 transformation temperature (heating) by fitting a power law relation 
to the experimental data reported in [27l [43] , see appendix A of [45] . 

As an example, for simulation of non-isothermal phase transformation, we have 
calculated the fraction of volume transformed during a j3 transition of Zircaloy-4 
material by solving equation ([4]) employing a similar triangular temperature history 
shown in figure [21 but extending its tail to 75 s. The initial condition utilized for 
heating is y(0) = 0, while on cooling y(tf) = 1, where tf corresponds to the end time 
of heating or the beginning of cooling. The Runge-Kutta method of order 4 and 5 
[46] was used to integrate equation ([4]) . Figure [4] shows y as a function of time and 
figure [5] depicts y versus temperature. In figure [5] we have included the equilibrium 
curve using equation (TT8]) and the corresponding experimental data reported in [27] . 
The computations on heating/cooling at the rates ±10 K/s are in agreement with the 
data reported in [57]. We note the asymmetry between the heating and the cooling 
in the computations. This is a reflection of experimental data from which the model 
has been adjusted to, which in turn emanates from the thermodynamic instability 
of the system. A more detailed comparison between experimental data and model 
computations, at several heating/cooling rates for zirconium alloys, is presented in 



4. Discussion 

The method utilized here to calculate the volume fraction of the favoured phase as 
a function of time and temperature can be applied to any input thermal history 
since we solve equation (0]) numerically. The question that may arise is how well 
this approximate equation compares with the prevailing exactly solved models in 
isothermal conditions. 
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Figure 4. Calculated volume fraction of /3-phase in Zircaloy-4 versus time during 
heating and subsequent cooling at the rates of ±10 K/s, using an extended version 
of the temperature history shown in figure [2] 




Figure 5. Calculated volume fraction of /3-phase in Zircaloy-4 versus temperature 
during heating and subsequent cooling at the rates of ±10 K/s, using an extended 
version of the temperature history shown in figure [2] The equilibrium curve is 
calculated using equation (1 1 8 I t and the corresponding measured data are from 
Forgeron et al. (2000). 
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Figure 6. Calculated volume fraction y of /3-phasc in Zircaloy-4 versus its rate 
dy/dt during heating and subsequent cooling at the thermal rates of ±10 K/s, 
using an extended version of the temperature history shown in figure [2] 



For isothermal conditions, equation (|6|) yields 

<j> = 1 - e -KT)(t-to) j (21) 

where to is the incubation time for the onset of phase transformation. This relation 
is a special case of the Avrami model [7J I25j related to site saturation transformation 
on grain surfaces under isothermal conditions, namely 

2S 

^ = l-exp[-— v(t-t )}, (22) 

where S/V is the surface area to volume ratio and v the interface velocity (growth 
rate) of the nucleus. Moreover, as has been shown by Cahn [25], equation ([22]) is 
theoretically exact for the systems in which the nucleation rate is sufficiently rapid 
such that site saturation occurs early in the reaction period. Indeed, he calculated 
that the time at which the site saturation occurs is t s = (I s w 2 ) -1 / 3 , where I s is the 
nucleation rate per unit area of the grain boundary, see |Appendix A| Our considered 
model is related to a special case of the KJMA theory where nucleation rate is high 
and site saturation occurs early in reaction, see | Appendix A[ 

A noteworthy observation in our computations, which is a reflection of 
experiments reported in |27j . is the asymmetry in phase transformation between 
heating and cooling carried out at the same thermal rating, cf. figures 2] and [5J To 
illustrate this asymmetry, we have plotted the phase portrait of our computations 
in figure [5] This is a manifestation of the different initial conditions and the 
different paths of the two processes. In our case, we had to extend the cooling 
tail of the temperature history (cf. figure [2J to 950 K to complete the phase 
transformation, see figure [5J A thermodynamic explanation of this phenomenon even 
at a phenomenological level is hard for us to ascertain here; our model is purely a 
kinetic one. 

The asymmetry is also true for the thermal rate dependence of the start 
temperature of phase transformation, cf. equations in appendix A of [45] . This 
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Figure 7. Calculated volume fraction of /3-phase in Zircaloy-4 versus temperature 
during heating and subsequent cooling at the rates of ±0.1 Ks , using an 
expanded version of the temperature history shown in figure [2] This figure should 
be compared with figure [5] 



type of hysteresis is ubiquitous in first-order phase transitions in solids. It is a 
manifestation of the solid state supercooling effect observed, e.g. in zirconium 
alloys, even at slow coolings (47l [48] . One reason for this behaviour is that the 
transformation process under non-isothermal conditions is away from equilibrium and 
is accompanied by dissipation of energy. The hysteresis should decrease with the 
decrease in transformation rate and should disappear at infinitely slow transitions. 
In other words, at infinitely low heating/cooling rates, the temperature variation of 
y should collapse on the equilibrium curve. To check this, we have repeated the 
aforementioned computations with a thermal rate two orders of magnitude lower, i.e. 
at ±0.1 Ks -1 by expanding the thermal history over 6000 s. The results are depicted 
in figure[7J which should be compared with figure^ We should, nevertheless, mention 
that in literature, the start temperature of first order phase transition has been related 
to the temperature dependence of the incubation time for nucleation [49j [45] . 

5. Summary and outlook 

In summary, we have studied the overall phase transformation kinetics in alloys under 
non-isothermal conditions using a differential analysis method. The method was 
applied successfully to the massive phase transformation in a zirconium alloy, for which 
the time/temperature variations of the volume fraction of the favoured phase under 
heating and cooling were calculated. A trait of non-isothermal or non-equilibrium 
phase transitions is the dependence of the onset of the transformation temperature on 
heating/cooling rates. We used best-fit empirical relations to account for this effect for 
Zircaloy. However, we alluded to the intimate relationship between this temperature 
and the incubation time for nucleation. This relationship warrants further analysis 
and may be studied in the context of transient nucleation phenomenon. 

Non-isothermal transformation kinetics is usually treated with the additivity rule, 
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where the transient heating/cooling is treated as a series of small isothermal steps. 
A formulation, which embraces both additive and non-additive situations, has been 
suggested. Moreover, a method based on the concept of path-integral, which accounts 
for all the possible thermal histories to reach the final state, has been presented. 
The path-integral approach has the potential to provide the optimal thermal path 
for phase transformation to attain the desired microstructure, and thereby material's 
macroscopic behavior. 
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Appendix A. Kolmogorov-Johnson-Mehl-Avrami model 

The overall phase transformation kinetics describes the time/temperature evolution 
of the volume fraction of the newly transformed phase onto the stable phase. The 
transformations are supposed to occur by nuclcation and growth mechanisms (i.e. 
first order transition). A mean field model was formulated and solved exactly by 
Kolmogorov [4], Johnson and Mehl [5], Avrami [6], Evans [50], and Jackson [51 j . 
apparently all independently. The basic assumptions of the model are : (i) the system 
is infinite in extent, and hence, boundary effects are neglected; (ii) nucleation is a 
stochastic process and on average takes place uniformly; and (iii) the growth of new 
phase grains ceases at the mutual points of their contact, but continues vigorously 
elsewhere; see ref. [52] for a recent appraisal. 

Kolmogorov [4j rigorously derived a relation for the time evolution of the 
transformed phase <f> in three dimensions. In a d-dimensional space, we write 

cb(t) = 1 - e~^\ (A.l) 

where <f> e (t) is the so-called extended volume fraction, expressed as 

Mt)=C d [ I(s)R d (t,s)ds, (A.2) 



Jo 

where I(t) is the nucleation rate (per unit volume and time), R(t,s) is the radius of 
the new grain (phase) at time t that was nucleated at time s, and Cd is a shape factor 
given for a hypersphere by 

27T d / 2 

Cd = dfm' (A - 3) 

with r(») denoting the usual gamma function. The function 4> e (t) is the sum of 
volumes of all the transformed grains divided by the total volume of the system, 
supposing that the grains never cease growing and the new ones keep nucleating at 
the same rate in the entire material. 

For R >> R* 7 where R* is the critical radius for nucleation, R(t,s) can be 
expressed in terms of the growth rate of the new phase as 

R(t,s) = [ v(r)dr, (A.4) 
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where v is the interfacial growth velocity. If now, the nucleation rate and the interfacial 
velocity assume anomalous power law dependencies on time, namely, I(t) ~ Iot a and 
v(t) ~ VQt 13 , then using equation (|A.4[) . equation (|A.2[) can be expressed as 

m = C (A.5) 



'3 

where m is an overall growth exponent of the favoured phase 

m = 1 + a + (1 + p)d, (A.6) 
and T g is a characteristic time for the growth of the favoured phase 

1 l/m 



(1 + /3) 



d+l 



c d i v$r(i + d)\i + f3Ji+d 



(A.7) 



with (z)„ = T(z + n)/r(z) being a Pochhammcr symbol. 

Here, we have shifted the time origin as t => t — fj„ c , according to modelling of 
transient nucleation in [53] . where U nc is called the incubation time and is found to 
be 

t $* 

Una = t* In-—, (A.8) 

where t* = R* /vq and $* is the energy barrier for nucleation in kelvin; see also 
Iwamatsu's calculations for d = 2 [54]. To make simplifications, recall the classical 
nucleation theory and assume that Iq is the steady-state nucleation rate with the form 
Iq = I* cxp(— <£>*/T), i.e. the number of critical nuclei per unit volume and time; e.g. 
see chapter 10 in [55]. Hence, according to equation (|A.7j) . r g = t* exp($*/?7i,T). 
Putting t ~ T g in equation (|A.8[) gives 

$* 

f — . (A.9) 
mi 



* 

ml 



In In — 

T T* 



So, as T is lowered, both r g and ti nc are increased, albeit their unlike temperature 
dependence. We note that in the aforementioned formalism with a = 0, (3 = and 
d = 3, with m = 4, the frequently used expression for the KJMA model is recovered. 

Cahn |25] derived relations for (f> e for particular cases where the new phase in a 
polycrystalline material nucleates on grain boundary surfaces, grain edges, or grain 
corners, under the condition that it grows with constant velocity. Two limiting cases 
corresponding to high and low nucleation rates (relative to growth rates) were found. 
At high temperatures, the nucleation rate is low and site saturation may not occur. 
If again / ~ I$t a and assuming random nucleation on either grain surfaces, edges, 
corners, or interiors, Cahn found (for d = 3) 

^W0k h4tm - (A ' 10) 
Moreover, since Iq would be proportional to the number of available nucleation sites, 
then Jo oc l/G d ~ J , where G is the parent phase grain size, j — d — 3 for corner, 
j = d — 2 for edge, j = d — 1 for surface, and j = d for intragranular nucleation. 
Equation (|A.10[) , except for a numerical factor and a shift in time origin, is close to 
the general equation (j A.5[) . 

For high nucleation rates, site saturation occurs early in the reaction and <f> e 
becomes independent of nucleation rate. The reaction is completed when t = 0.5G/v 

4> e - C nPn (vt) n . (A. 11) 
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Here, n = 1, 2, 3 correspond respectively to site saturation on grain surfaces, 
edges and corners; and p n denotes the ratio of the respective surface area, grain- 
edge length and grain corner number per unit volume; note that 4> e is independent 
of nucleation rate and is driven by the growth rate. If we put v(t) ~ vot@ , then 
e ~ C nPn v%t< 1+ ft = M" (1+/3) . 

Thus, for n = 1 (grain boundary site saturation) and [3 — equation (jA.lj) can 
be written as 

4> = l-er kl \ (A.12) 

where kx = CipiVa = 2(S/V)vo, which except a shift in the time origin, is exactly 
the same as equations (|2i p -(|22 |) , and is equivalent to the model we have used in our 
analysis. 

It is worthwhile to mention that Ham |56j obtained a similar kind of relation as 
(|A.12j) for diffusion-limited precipitation of second phase in a crystalline. In Ham's 
model, one may imagine that a polycrystalline consists of an array of parallel cylinders, 
in which the solute atoms precipitate on cylinders' surfaces. For this model, the 
precipitated fraction of excess solute w can be calculated as a function of time t fairly 
accurately by the simple relationship 

w = l-e"* /T0 , (A.13) 

where t = £ 2 /(alD), I is the inter-precipitate distance, a the lowest eigenvalue of 
the prevailing diffusion equation, and D is the solute diffusivity. Hence, the theoretical 
basis of the model used for calculations of phase transformation kinetics in this paper 
may be attributed to a special case of the KJMA formulation, or diffusion-limited 
phase transformation on surfaces. 
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